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We investigate whether the Affleck-Dine (AD) mechanism works when the contribution of the two- 
loop thermal correction to the potential is negative in the gauge-mediated supersymmetry breaking 
models. The AD field is trapped far away from the origin by the negative thermal correction for 
a long time until the temperature of the universe becomes low enough. The most striking feature 
is that the Hubble parameter becomes much smaller than the mass scale of the radial component 
of the AD field, during the trap. Then, the amplitude of the AD field decreases so slowly that 
the baryon number is not fixed even after the onset of radial oscillation. The resultant baryon 
' asymmetry crucially depends on whether the Hubble parameter, H , is larger than the mass scale 

I of the phase component of the AD field. Me, at the beginning of oscillation. If _ff < Me holds, the 

. formation of Q balls plays an essential role to determine the baryon number, which is found to be 

^ ' washed out due to the nonlinear dynamics of Q-ball formation. On the other hand, ii H > Me 

^ holds, it is found that the dynamics of Q-ball formation does not affect the baryon asymmetry, and 

, that it is possible to generate the right amount of the baryon asymmetry. 

Q ■ PACS numbers: 98.80.Cq 

cn \ 

CN ' I. INTRODUCTION 

> : 

1^ \ Baryogenesis is one of the main issues in the theories of the early universe. There have been proposed various 
I mechanisms, among which AfHeck-Dine (AD) mechanism is a promising candidate in theories with supersymmetry 
P\| (SUSY). One of the features which distinguish supersymmetric theories from ordinary ones is the existence of fiat 
directions, of which the AD mechanism makes use. In fact, there are many fiat directions in the minimal supersym- 
metric standard model (MSSM), along which there are no classical potentials. Such a flat direction is described by 
, a complex scalar field <i>, called AD field. It is crucial to determine the dynamics of the AD field to estimate the 
resultant baryon asymmetry. However, it has been found that the AD baryogenesis is complicated by the two causes. 
Q-i One is the thermal effect |3, 0, : and the other is the existence of Q balls ^^J^. 

As for the thermal effects, we must take account of those which appear at both one-loop and two-loop orders. The 
^ former contribution to the effective potential of the AD field is the thermal mass term which might cause the early 
oscillation of the AD field. The latter effect, on which we would like to focus in this paper, induces the logarithmic 
^ ' potential (which we call thermal log potential hereafter) and hence totally changes the dynamics if its coefficient is 
• ^ , negative. In that case, the AD field might be trapped far away from the origin for a long time until the temperature of 
' the universe becomes low enough. In particular, such a trap is likely to occur in the gauge-mediated supersymmetry 
H [ breaking (GMSB) models 's'l where the potential for the AD field becomes flat at large amplitude. Therefore, we 
assume the GMSB models in this paper. Due to the enduring trap by the negative thermal log potential, the dynamics 
of the system is so involved that little has been known about whether the AD baryogenesis really works. Although 
the sign of the thermal contribution to the potential at two-loop order depends on the choice of flat directions, it is 
expected to be negative for many flat directions. Hence it is important to consider the case when the negative thermal 
log potential is effective. 

The other element which we need to incorporate is the existence of Q balls 0, 0, Q ■ After the AD field starts to 
oscillate, it experiences spatial instabilities and deforms into nontopological solitons, Q balls. Especially, the Q-ball 
formation is inevitable for a flat potential as in the GMSB models 0,0,^3- Once Q balls are formed, they absorb 
almost all the baryon numbers |^llCll|. In the usual picture, the Q-ball formation takes place after the baryon number 
is fixed. However, as shown later, the Q-ball formation plays an essential role to fix the baryon number if the AD 
field is trapped by the negative thermal log potential. 

In this paper, therefore, we consider the AD baryogenesis with the negative thermal log potential in the GMSB 
models, taking into account the Q-ball formation. Since the negative thermal log potential retards the oscillation of 
the AD field considerably, the dynamics is quite different from that for the ordinary AD baryogenesis, and we need 
to examine the scenario carefully. It is found that the existence of the Hubble-induced A-term spoils the baryogenesis 
through the Q-ball formation, while it is possible to explain the baryon asymmetry of the present universe in the 
absence of the Hubble-induced A-term. 
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II. AFFLECK-DINE MECHANISM AND Q-BALL FORMATION 



In this section we briefly review the Affleck-Dine mechanism and properties of the Q baUs. In the MSSM, there 
exist flat directions, along which there are no classical potentials in the supersymmetric limit. Since flat directions 
consist of squarks and/or sleptons, they carry baryon and/or lepton numbers, and can be identified as the AD field. 
These fiat directions are lifted by both the supersymmetry breaking effects and the nonrenormalizable operators with 
some cutoff scale. In the gauge-mediated SUSY breaking model, the potential of a flat direction is parabolic at the 
origin, and almost ffat beyond the messenger scale 0, fill ITll| , 

m « Ms), 

Vgauge^i ,,4 f, n^,.. .. . (1) 

M2 



M^(iogif^) m->Ms), 

s ■ 



where is a soft breaking mass ~ 0(1 TeV), Mp the SUSY breaking scale, and Ms the messenger mass scale. Since 
the gravity always exists, flat directions are also lifted by the gravity- mediated SUSY breaking effects [T^ . 



^grav — ^3/2 



1 -I- iv: log 



1^ 



(2) 



where K is the numerical coefficient of the one-loop corrections and Mq is the gravitational scale (~ 2.4 x 10^^ GeV). 
This term can be dominant only at high energy scales because of small gravitino mass ^ 0(1 GeV). 

The nonrenormalizable operators, in addition, lift fiat directions. Assuming that there exist all the nonrenormaliz- 
able operators consistent with both the gauge symmetries of the standard model and the R-parity in the superpotential, 
most of fiat directions are lifted by the following operator: 

W = ,1 , (3) 

where M is the cutoff scale. Not only does it lift the fiat potential as 

= ' (4) 



but it supplies the baryon number violating A-term as 



^^ = «™;^*" + H-C., (5) 

where Um is a complex constant of order unity, and we assume the vanishing cosmological constant. The schema of 
the zero-temperature potential Vgauge + Vgrav + Vnr. is shown in Fig. ^ 

For the AD baryogenesis to work successfully, it is necessary that the AD field has a large expectation value during 
the infiationary epoch [23|. To this end we require the existence of the four-point coupling of the AD field to the 
inflaton / in the Kahler potential, which leads to the negative Hubble-induced mass term: 

= -c/fi72|$|2, (6) 

where cr is a positive constant of order unity. Also there might be a Hubble-induced A-term, which comes from the 
three-point coupling of the AD field to the inflaton in the Kahler potential. We write it as 

yAH=aH^^,'^- + ^.c., (7) 

where an is a complex constant of order unity. Though the presence of this term usually does not change the resultant 
baryon asymmetry, it does change the dynamics of the AD field when the negative thermal log potential traps the 
AD field, as shown later. Since the Hubble-induced A-term is not necessarily present, we will consider both cases 
with and without this term. 

Let us consider the thermal effects on the effective potential for the AD field. During the oscillation of the inflaton, 
there is dilute plasma with temperature T ~ {T'^jjHMg)^^^, where Trh is the reheating temperature. If the AD 
field directly couples with fields ipk in the thermal bath, it acquires a thermal mass term in the effective potential at 
one-loop order: 

/fc|*l<T 



/fi)T2|ci>|2 for |$| < /(-It, (8) 
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where Ck is a real positive constant of order unity, and fk is a Yukawa, or gauge coupling constant between the AD 
field and tpk H,^!- For simplicity, we omit the subscript k hereafter with the understanding that /(i) denotes either 
Yukawa or gauge coupling constant. Note that this effect is exponentially suppressed by the Boltzmann factor when 
the effective mass of the thermal particle, is larger than the temperature. In addition to this term, there 

is another thermal effect on the potential, which appears at two-loop order, as pointed out in Rcf. J^. This comes 
from the fact that the running of the coupling constant is modified by integrating out heavy particles which directly 
couples with the AD field. This contribution is given by 0| 

T4^)~aT/(%TMog^, (9) 

where qt is an 0(1) constant, and f^2) = g^iV^iU'^g'^ {d and y are the gauge and Yukawa coupling constants). The 

sign of Or depends on the flat directions, and it is expected to be negative for many flat directions, for instance, udd 
direction, though it is complicated to determine the precise value of a^. Without specifying flat directions, we would 
like to focus on the case with the negative ot in the following discussion. (The AD baryogenesis for the positive ar 
was considered in Re f. 10].) 

When the AD fleld starts coherent rotation in the potential, the baryon number is generated and fixed immediately 
in the ordinary AD mechanism, since the A-terms becomes ineffective soon after the oscillation of the AD field. Then 
the baryon number density is estimated as 

nsitosc) ^ euj(t)l^^, (10) 

where e(^ 1) is the ellipticity parameter, which represents the strongness of the A-term, and uj and iposc are the 
angular velocity and the amplitude of the AD field respectively, at the beginning of the oscillation (rotation) in its 
effective potential. However, this is not the case when the AD field is trapped by the negative thermal log potential 
given by Eq. 0, which will be discussed in the next section. In fact, the baryon number is not fixed even after 
the AD field starts coherent oscillation, if the Hubble parameter is much smaller than the mass scale of the radial 
direction. It should be emphasized that Eq. (|10|) is valid as long as the baryon number is fixed soon after the oscillation 
commenced. 

Lastly, we comment on the Q balls. In the course of the coherent oscillation, the AD field experiences spatial 
instabilities, and deforms into nontopological solitons, Q balls When the effective potential is dominated by 

the zero-temperature term Vgauge , the gauge- mediation type Q balls are formed, whose properties are as follows : 

Mq ~ MfQ^/\ Rq - Mp^Q^/\ (11) 

where Mq and Rq are the mass and the size of the Q ball, respectively. If the mass per unit charge, MfQ~^/'^, is 
smaller than the proton mass ~ IGeV, the Q ball is stable against the decays into nucleons, from which it follows 
that Q balls with very large Q can be stable. From numerical calculations 0, 0| , it is known that Q balls absorb 
almost all the baryon charges which the AD fleld obtains, and the typical charge is estimated as ^3 



^ 4 

^osc 



where /3 « 6 x 10^*. Consequently, the present baryon asymmetry should be explained by the charges which come 
out of the Q balls through the evaporation and decay of Q balls. 

In the case of the unstable Q balls, they decay into nucleons. The decay rate is given by 0| 



dQ 



dt 



<^ (13) 



192^2 ' 



where A is a surface area of the Q ball. 

In the case of the stable Q balls, the evaporation is the only way to extract the baryon charges from Q balls. The 
total evaporated charge from the Q ball is estimated as [Tfl llS > 

AQ ^ 10^5 (^) ( ) gVi2. (14) 



Hence the baryon number density is suppressed by the factor AQ/Q, in comparison with the case of unstable Q-ball 
production. 
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III. DYNAMICS OF AFFLECK-DINE FIELD 



In this section, we present a detailed discussion on the dynamics of the AD field. First we derive the condition that 
the AD field is trapped by the negative thermal log potential. In the following argument we set \am\ = ch — |a_ff I — 
Cfc = —aj- = 1 for simplicity. During inflation the Hubble parameter is roughly constant, and the initial amplitude of 
AD field is set far away from the origin due to the negative Hubble-induced mass term Eq. The position of the 
global minimum is 

~ (iJM"-3)i/("-2) ^ (15) 

where (j) = •\/2|*I'|. After inflation ends, the Hubble parameter decreases, and so does the minimum of the potential. 
Since the AD field tracks the instantaneous minimum, Ea. (|15(l still holds after inflation. The AD field is trapped by 

the negative thermal log potential when |Vy | > \Vh\ is satisfied. Here we assumed that the trap due to the negative 
thermal log potential occurs while the AD field tracks the instantaneous minimum given as Ea. (|15|l . The temperature 
Tc and the amplitude (j)c at the beginning of the trap are 



f^ - T^^ Mg- M-— for Trh < n 

for Trh > T* 



a,T^M-^)-\ (16) 



where we defined 



T^ = !^-^H'I^' M-^. (17) 



It should be noted that Trh > requires a very high reheating temperature which might be constrained by the 
gravitino problem |lll |. However, if the Q balls dominate the universe, they must decay and reheat the universe with 
the decay temperature, Td i$ O(lGeV). Thus the gravitino problem is alleviated due to the entropy production of the 
decay of Q balls in this case, and such a high reheating temperature might be possible. If it were not for the gravitino 
problem, the reheating temperature could be as large as Tinf = 10^^ GeV. We take this value, because the COBE 

data implies that vlJ^/eli'^ ~ 6.7 x lO^^ GeV 13, where e^r < 1 is the slow-roll parameter, and the instantaneous 
reheating is assumed for conservative discussion. 

Next we consider when the trap due to the negative thermal log potential ends. The AD field starts an oscillation 
when \Vj, \ < Vgauge is first satisfied, hence the temperature Tosc and the amplitude (f>osc at that moment are 

Tosc = f{2)^^P ' 

(j)osc = min [cj)a, 0/3] , (18) 

1 

where = (MJ^M"^"^) ""^ and (t>p = M|,/to3/2. Note that the amplitude of the AD field, (/)osc, at the beginning of 
oscillation is different from that in the usual scenario. In the ordinary AD mechanism, the AD field starts oscillating 
at 4>a or (j)^ = (to3/2M"^^)'^^ depending on the form of the potential shown in Fig. ^ if one does not take account 

(2) 

of the thermal effects. Such modification is due to the different dependence of Vh and Ij. on 

From the above discussion we can derive the condition that the AD field is trapped by the negative thermal log 
potential. It requires the following inequalities to be satisfied. 

0c > Max [0a, 07], (19) 
/(i)0c > Te, (20) 



> Ms. (21) 



Then we arrive at the following conditions. 



10"^ GeV< m3/2 <m'fj>^, 
10TeV< Mf <M^p\ 

Tli^H< Trh <tII'^, (22) 
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where rn^^j^.^^ ^f^\ and are given as 



(«) ^ 

^3/2 — 



M 



(u) _ 




for Trh < T^, 



1 GeV, f^^-'^M^-^M-^''-^^ 



-2 2(ti-4) 



(23) 



for Tj^ff > T„ M > /(i^V(2r' ^G, 

for T/jH < 



(24) 



for Tj^H > T„ A/ > / "fV(lf '"^G, 



^RH — 




for M>/(f'/(,) 



-2 2(71-4) 
-3 11-3 



for M 



-2 2(ti-4) 
-3 £ 71-3 



(25) 



(26) 



M, 



G, 



where we do not take account of the gravitino problem, and the bounds for Mp and 7713/2 include the constraints 

due to the scheme of the GMSB. Notice that T^p^ comes from > Max[(/)c, 0^] and T^^^ is obtained from the 
condition /(i)0c > ^c- Also, we must impose the following condition in order to take the gravitino problem into 
consideration. 



lOOGeV - ITeV 



for 2/i^kcV< 



< lOOkeV 



Trh < 7 X 



lOTeV xh^ ( "^^^^ ) ( ^ for TO3/2 > lOOkeV, 
VlOOkeV/ VlTeV/ ^' 



(27) 



where h is the present Hubble parameter in units of lOOkm/sec/Mpc, togs is the gaugino mass for SU{3)c, and 
7 reflects the entropy dilution due to the decay of Q-balls. If Q balls are subdominant component of the universe 
throughout, we have 7 = 1. On the other hand, if Q balls dominate the universe, it is given as 



Trh 



for Te„>T, 



RH^ 



for T,, < T, 



(28) 



RH, 



where Teq is the temperature at which the energy of Q balls becomes comparable to the total energy of the universe, 
and Td is the decay temperature of Q balls given as 



90 



IMfMg 
A8nQi 



(29) 



where 5, counts the effective degrees of freedom of the radiation. 

So far we have outlined the dynamics of the amplitude of the AD field. Now we would like to focus on the motion 
of its phase, 9 = arg[<i>], which essentially determines the baryon asymmetry. The question we have to ask here is 
when the baryon number is fixed. For this purpose let us begin our analysis by following the evolution of the baryon 
number density. For a moment, we do not consider the effect of the Q-ball formation. The equations of motion for 
and are 



• dV 

3H<j) + TTT - 



= 0, 



3H + 2-\e 



(30) 
(31) 
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The baryon number density is given as 

= qO<P\ (32) 
where we defined that the baryon charge of the AD field is q. Its evolution is governed by the following equation, 

dV 

UB + iHuB + = 0. (33) 

The first and second terms mean that the number density is diluted by the expansion of the universe, and it is the third 
term that represents the generation of the baryon number. Note that the only A-term, VA(+Vyi//), can contribute to 
the third term. 

In the first place we review how the baryon number is fixed in the ordinary AD mechanism. In this case, the 
AD field starts the oscillation when the Hubble parameter becomes equal to the mass scale of the radial direction, 
M^. Thereafter the amplitude is reduced by the Hubble frictional term. Accordingly the magnitude of the A-term 
decreases rapidly, since it is a positive power of </). Then it can be shown that the A-term cannot generate the 
baryon number efficiently after the AD field starts the oscillation as follows. While the AD field tracks the slowly 
decreasing minimum given as Eq. (|15l) . the motion of 6 is strongly damped by the frictional term because (t)/(f) ^ H 
and H ~ ^ Me, where Mg is the mass scale of 6. The second equality is satisfied if the Hubble-induced A-term 
exists. Hence the characteristic time scale for 6 is (Alg / H)''^ , and the baryon number density is expressed as 

nB=qH'-q^^'. (34) 

Then the evolution equation for is given as 

dV 

TiB + SHuB = ~q^, 
oO 

~ Hub, (35) 

where we used Eq. ()34|l in the last equality. From Ea. l|35|) . we can see that the baryon number continues to be 
critically generated while the AD field is trapped by the negative Hubble mass term. That is to say, the decrease 
of the baryon number density due to the expansion of the universe is compensated by the baryon number violating 
term to some extent, as long as the AD field tracks the instantaneous minimum given as Eq. (|15|l . However, once 
it starts the oscillation, the magnitude of the A-term is significantly reduced, and also the time scale for 9 becomes 
much smaller : '-^ ^ H~^. Thus the generation of the baryon number stops soon after the AD field starts 
oscillating in the ordinary AD mechanism. 

The usual picture outlined above does not apply to the case when the AD field is trapped by the negative thermal 
log potential. First we assume the existence of the Hubble-induced A-term, and see how it spoils the baryogenesis. 
The relations among the mass scales, M^, Mg and H, are then Af^ ~ Mq ~ H before the trapping due to the thermal 
effect. Once the AD field is trapped by the negative thermal log potential, its amplitude, ^ (/j^2)^^-^"~^)^^'"~^''> 
decreases relatively slowly. While both and Mg are decreasing functions of slowly varying 0, the Hubble parameter 
diminishes faster. Thus, M^ ^ Alg '3> H holds until the trap ends. Since the mass scale of 9 is larger than the Hubble 
parameter, 9 starts oscillating around its minimum which is determined by the superposition of Va and Vah, as soon 
as the AD field is trapped. During the trap, the amplitude of the oscillation of 9 decreases continuously |2j, so that 
9 stays very close to the bottom of the valley of the A-term when starts an oscillation. Therefore the AD field 
rotates the orbit with very large oblateness, as shown schematically in Fig. |21 It should be noted that the amplitude 
of the oscillation of (p remains almost constant due to the small Hubble parameter, which has a decisive influence on 
the evolution of the baryon number density. In order to follow the evolution of ub, we estimate the third term in 
Eq. I|33|) . which is now given as 

dV 

~ qM^^ = 9-^MSnB, 

MgUB , (36) 



< 



where the last equality holds at the end point of each oscillation. Therefore the baryon number density continues to 
be generated around the end point, since Mg ^ H. Put simply, the valleys of the A-term become effectively deeper 
due to the enduring trap, which causes the intermittent generation of the baryon number. 
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Consider now the implication of the Q-ball formation. The AD field experiences instabilities when it oscillates 
in the Vgauge dominated potential, and deforms into the gauge- mediation type Q balls. Since the growth rate of 
the instability is of the order of ^ M^/cjiosc Q balls are formed soon after the AD field starts a radial 
oscillation '2^. The most important part of this argument is that the baryon number is not fixed during the process 
toward the Q-ball formation. When the fluctuation becomes comparable to the homogeneous mode, the system goes 
into the nonlinear stage, and the AD field is kicked back and forth at random in the phase direction by the A-term. 
It is not until Q balls are formed that the baryon number is fixed. The reason is as follows. It is known that Q balls 
absorb almost all the baryon number. Inside Q balls, the AD field rotates a circular orbit, so the average of the third 
term in Eq. (|33() over a round of the orbit vanishes: 

dV 

which means the baryon number is fixed. Since the baryon number evolves indiscriminately until the Q-ball formation 
is completed, the resultant baryon asymmetry is affected by the highly nonlinear physics of the Q-ball formation. In 
rough estimate, the created baryon number before the system goes into the nonlinear stage is 



n^)~-^M,VLsinn0osc, (38) 



where 0osc(^ 1) is the phase amplitude at the start of the radial oscillation, and ciM^ ^ with ci ~ log [(ji/Scj)) = 0(10) 
is the time for the fluctuation to become nonlinear. When the fluctuation of the phase becomes large, i.e., 59/6 ~ 0(1), 
the AD field moves beyond the potential hill in the phase direction and the baryon number production restarts with 
random phase, 9i{x.) ^ 0(1). Here 9i depends on the position x. Hence this process induces the fluctuation of the 
baryon number density, and it takes place over and over again until Q-ball formation is completed. If the process 
lasted eternally, however large the baryon asymmetry is, it would be completely erased due to the baryon number 
violating process in equilibrium. Note that the baryon number is odd under the CPT transformation. In fact, this 
process is effective for a finite time, so that the only finite baryon asymmetry is annihilated. We can estimate the 
magnitude of the baryon asymmetry, which is erased due to the Q-ball formation processes, by calculating the typical 
value of the fluctuation of the baryon number density : 



1 

SriBi^) — M^(j>l^^'^ sin n9^{x), 



~ ±^^e<^L, (39) 

where C2 ^ 0(10) is the number of the restarts, and the contribution of the random phase 9i is evaluated as random 

walk with C2 steps in the second line. Since 9osc ^ 1 in the present situation, we have n^' <^ Sub, which means that 
TT-B is dominated by the fluctuation due to the Q-ball formation process. Moreover, the dynamics is chaotic in the 
sense that a tiny change results in a quite different baryon number. Thus, the baryon asymmetry generated during 
the homogeneous evolution is swallowed up in the large fluctuation, leading to the baryon symmetric universe. 

We performed the numerical simulation in order to confirm whether the baryon asymmetry is really determined 
through the Q-ball formation. For the numerical calculation, we take the variables dimensionless as follows. 



Ms 

X. = (40) 



The initial conditions of the homogeneous parts are taken as 

Re[v3(0)] = 1.0 X 10^ 
Im[v3(0)] = 1.0 X IQ-^ 



ip'{Q) = 0.0, (41) 



where the prime denotes the derivative with respect to r, and we redefined the phase of the AD field so that the 
real direction coincides with one of the valleys of the A-term with the mass scale Mg/M^ = lO""'^'^. Taking account 



8 



of fluctuations which originate from the quantum fluctuations, we give ip fluctuations with amphtudes 10~ times 
smaller than the homogeneous mode. We have confirmed that the smaller fluctuations just delay the formation of Q 
balls. For simplicity we neglected the cosmic expansion. The results are shown in Fig. 13 It shows that the baryon 
number is generated at the end point of each oscillation for some time ~ ciM^^, but it soon goes into the chaotic 
phase until the Q-ball formation is completed. If we change the initial seed for the fluctuation, it is probable that 
the resultant baryon asymmetry inside the box of the lattices becomes positive. The reason why the resultant baryon 
number takes nonzero value is that the size of the box under consideration is finite. If the same calculation were done 
with the infinitely large box, the resultant baryon number would be exactly zero. 

Lastly we comment on the case without the Hubble-induced A-term. Now > H > Mg at the onset of the 
oscillation can be realized since the mass scale of 9 is smaller than that in the previous case. Then the initial phase 
amplitude 9osc is comparable to the random phase 9i, hence the net baryon number can be nonzero. In the next 
section wc estimate the resultant baryon asymmetry in this case, taking account of the Q-ball formation. 



IV. BARYOGENESIS 



As we have shown in the previous section, the Hubble-induced A-term spoils the baryogenesis through the Q- 
ball formation. On the other hand, if it were not for the Hubble-induced A-term, we can avoid this problem since 
Mcf, > H > Mg at the beginning of the oscillation is possible. Then the baryon number just before the Q-ball 
formation is completed is roughly estimated as 



(42) 



where ci ^ C2 ^ 0(10), and we used 9osc ^ 0{1) in the present case. Thus we expect that the finite baryon asymmetry 
remains even if the AD field gets randomly kicked in the phase direction during the nonlinear stage. In other words, 
the resultant baryon number is almost determined due to the dynamics before the AD field goes into the nonlinear 
phase, so subsequent Q-ball formation does not affect the baryon asymmetry. 

To satisfy the requirement that the net baryon number is generated, the following condition must be satisfied: 



Mg < Ho 



(43) 



where 



Mg 



Ho 



Then the baryon number is generated as 



' 'rn^/24'osc 



T^hMg 



I Mc 



nsitosc) 



for Tosc > Trh, 



for Tosc < Trh- 



(44) 



(45) 



1 "^3/2 Cs, 

M^ M"-3 



(46) 



where it should be noted that = Mp/(j)osc is used instead of the Hubble parameter, because the kick due to 
the A-term lasts for St ^ M^^, and we also set the numerical factor to be unity. The almost all baryon number is 
absorbed into the subsequently formed Q balls with the typical charge given as Eq. (|12|l 23]. Hereafter we consider 
the two cases, depending on whether the Q balls dominate the universe or not. If the Q balls dominates the universe, 
they must be unstable and decay before the epoch of big bang nucleosynthesis (BBN). Thus, the mass per unit charge 
and the decay temperature are constrained as 



T. 



90 



MfQ- 



IMfMg 



> IGeV, 



> lOMeV. 



(47) 
(48) 
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Since the universe is reheated again, the gravitino problem is alleviated in this case. Accordingly the constraints 
Eqs.|(22J) and H27|) with 7 given as Eq. (|28|) are adopted. The baryon-to-entropy ratio is given as 

riB _ 3 nB{tosc)Td 
s ^ 4 M| 

30 X (|^)"*m37+^A4^"-''M-("-3)M| for m^^ > (Ml^-H.r^"-^^ 



(49) 



30x(|^) ~^ mz,2Mp--'A'r^M^ for 7713/2 < (a/^^-^M- ("-3) 



The allowed regions for 7713/2, Mp and Tjih satisfying the above constraints are shown in Fig^landEl where we set 
/(I) ~ /(2) ~ 0.3, M = Mg and n = 6. It shows that the baryon-to-entropy ratio can be large enough to explain the 
present asymmetry. For instance, it takes the value as 

-=^^^^ (10) ( lOOM^j U X lO^GeV j [m^) ' ^'^^ 

where 7773/2 and Mp is assumed to satisfy all the constraints stated above, so their values must be included in the 
allowed regions shown in Fig^and|Sl Though the above value is a little larger compared to the expected one ~ 10"^", 
it might be suppressed by the unknown CP phase ^oscwhich we assume to be of the order of unity here. Although 
Af > Mq is possible since the definition of M in Eq. Q includes a possibly small coupling constant, it does not work 
to take larger M. This is because the requirement that Q balls decay before the BBN sets the lower bound for Mp 
which is a positive power of M, so the allowed regions shown in Fig^ and |S1 would disappear for M ^ Mq- 

Next we consider the case that the energy of Q balls is subdominant before the BBN. If Q balls are unstable, they 
must decay before the BBN, otherwise their decay product might destroy the light elements and spoil the success of 
the BBN. If Q balls are stable, the resultant baryon asymmetry is suppressed by AQ/Q, where AQ is the evaporated 
charge. Also the abundance of remnant Q balls are constrained in order not to overdose the universe. Namely, 

MpQ^Qb Q 
77i„ e AQ 

where 777„ is the nucleon mass, and f^Q(B) is the ratio of the energy density of Q balls (or baryon) to the critical 
density. The baryon-to-entropy ratio is expressed as 

UB / 45 \ HBitosc) 

' ' X Ki X K2, (52) 




where 

for Tosc < Trh, 

= <( (53) 

for Tosc > Trh, 

for unstable Q balls, 
K2 = { (54) 
AQ/Q for stable Q balls. 

Since the reheating temperature is now constrained by the gravitino problem, the constraints Eas. (|22|l and (|27|) with 
7 = 1 are adopted. The allowed regions for 7773/2, Mp and Tj^h satisfying the above constraints are shown in Fig|^ 
and[3 where we set /(i) = /(2) = 0.3, M = Mg and n — 6. It shows that the baryon-to-entropy ratio can be as large 
as ~ 10~^^, which is a little smaller than the favored value ~ 10"^". However, it can be cured by taking smaller 
Af, the GUT scale for instance. FiglHl and El represent the allowed regions for 7773/2, Mp and Trh in the case of 
/(I) — /(2) = 0.3, M = 10^^ GeV and 77 — 6. It shows that the baryon-to-entropy ratio can be now as large as 10^^° 
for unstable Q balls, although it is too small to explain the present asymmetry for stable Q balls. For instance, it is 
given as 



riB 
s 



3 X 10 



-11 (l^y ( ^3/2 X ( Mp / Tnn \ ( M y 

V2OO/ VlOOMeV/ ^2 X lO^GeVy \lQ^Ge\f ) yW^^GeY ) ' ^ ' 



for unstable Q balls. 
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V. CONCLUSION 

We have studied whether the AD baryogenesis works in the gauge-mediated SUSY breaking models, if the AD field 
is trapped by the negative thermal log potential which is known to be ubiquitous in the MSSM. The most striking 
feature is that the Hubble parameter becomes much smaller than the mass scale of the radial component of the 
AD field, during the trap. Then, the amplitude of the AD field decreases so slowly that the baryon number is not 
fixed even after the onset of radial oscillation. A further important point is that the AD field experiences spatial 
instabilities and deforms into Q balls. Therefore the baryon asymmetry can be washed out due to the random kick 
in the phase direction during the process of Q-ball formation. Indeed, it has been shown that the net baryon number 
vanishes if the Hubble parameter is smaller than the mass scale of phase direction, since the randomly generated 
baryon number surpasses the baryon number just before the nonlinear stage begins. This is always the case if the 
Hubble-induced A-term exists. Also we confirmed the chaotic behavior of the baryon number with use of numerical 
calculation. However, the Hubble-induced A-term does not necessarily exist. In fact, we have found that it is possible 
to generate the right amount of the baryon asymmetry in the absence of the Hubble-induced A-term, if the Q balls 
are unstable. Thus, it should be emphasized that the AD baryogenesis is still a viable candidate, even if both the 
thermal effect and Q-ball formation are fully taken into account. 
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FIG. 1: Schema for the zero-temperature potential Vgauge + Vgra + Vnr- If »7i3/2 < Mp ^ M "-i is satisfied, the potential 
takes the form shown in the upper figure. The Vgauge-dominant region connects to the non-renormalizable term dominant 

region at (f)a ~ (mI-M""^) . On the other hand, if 7713/2 > Mp^^ M~ "-1 is satisfied, the potential takes the form shown 

in the lower figure. The Vgauge-dominant region connects to the Igrau-dominant region at (fip ~ , and then it connects to 

the non-renormalizable term dominant region at 0^ ~ {m:j/2M" ) "-^ . 
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FIG. 2: This figure illustrates how the AD field behaves after the beginning of the oscillation, when it is trapped by the 
negative thermal log potential. The thick solid line is the orbit of the AD field. The dashed lines represent the hills of the 
A-term, and the AD field oscillates in the valley surrounded by them. At each oscillation, it is kicked into the direction shown 
as the broad arrow. 
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FIG. 3: The evolution of the baryon number during Q-ball formation is shown. As shown in the lower figure, it can be classified 
into three parts. First the AD field is kicked around the endpoint at each oscillation (see Fig. leading to the intermittent 
baryogenesis (interval I). Then the behavior of the baryon number becomes chaotic, as the fluctuations grows (interval II). 
Finally, it is fixed after the Q-ball formation is completed. The upper figure is the blowup of the evolution for the baryon 
number during the interval I. 
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FIG. 4: The allowed regions for m3/2, Mp and Trh in the case that the Q balls dominate the universe are shown. Here we 
set /{!) = /(2) = 0.3, M = Mg and n — 6. 
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FIG. 5: The baryon-to-entropy ratio is plotted for the allowed regions in Fig.Q] 
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FIG. 6: The allowed regions for 01,3/2, Mp and Trh in the case that the Q balls do not dominate the universe are shown. 
Here we set /(i) = /(2) = 0.3, M = Mq and n = 6. The open cicle and closed triangle denote stable and unstable Q balls 
respectively. 
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FIG. 7: The baryon-to-entropy ratio is plotted for the allowed regions in Fig. El The open cicle and closed triangle denote 
stable and unstable Q balls respectively. 
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FIG. 8: The allowed regions for 7713/2, Mf and Trh in the case that the Q balls do not dominate the universe are shown. 
Here we set /(i) = /(2) = 0.3, M = 10^^ GeV and n = 6. The open cicle and closed triangle denote stable and unstable Q balls 
respectively. 
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FIG. 9: The baryon-to-entropy ratio is plotted for the allowed regions in Fig. |H| The open cicle and closed triangle denote 
stable and unstable Q balls respectively. 



